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ABSTRACT 

We show that discs accreting onto the magnetosphere of a rotating star can end up in a trapped 
state, in which the inner edge of the disc stays near the corotation radius, even at low and vary- 
ing accretion rates. The accretion in these trapped states can be steady or cyclic; we explore 
these states over wide range of parameter space. We find two distinct regions of instability, 
one related to the buildup and release of mass in the disk outside corotation, the other to mass 
storage within the transition region near corotation. With a set of calculations over long time 
scales we show how trapped states evolve from both nonaccreting and fully accreting initial 
conditions, and also calculate the effects of cyclic accretion on the spin evolution of the star. 
Observations of cycles such as found here would provide important clues on the physics of 
magnetospheric accretion. Recent observations of cyclic and other unusual variability in T 
Tauri stars (EXors) and X-ray binaries are discussed in this context. 
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1 INTRODUCTION 

Many accreting stars show evidence of the effects of a strong stellar 
magnetic field regulating the accretion flow onto the star. Accret- 
ing X-ray pulsars, for example, show flux modulation on timescales 
between ~ 10~ 3 — 10 3 seconds, which is attributed to the mag- 
netic pole sweeping through our line of sight, illuminate d by mat- 
ter accreting along field lines dDavidson & Ostrikerl 19731) . In some 
pulsars, this probe of the star's period has shown an evol ution in 
the spin period of the star in time (e.g. lBildsten et al]|l997h . which 
is attributed to the interaction with the surrounding material. On 
the other end of the energy scale, T Tauri stars also often have 
strong magnetic fields (up to 1-2 kG), and show som e evidence that 
their period is regulated by the disc-field interaction ("Germ an" et alj 
l2008h . 

From the earliest studies of stellar accretion researchers rec- 
ognized that a strong magnetic field c ould substantiall y influ ence 
an accretion disc surrounding the star. IPringle & Reesl dl972[) es- 
timated the 'magnetospheric radius' of the disc - the radius at 
which the magnetic field would truncate the disc - by equating 
the ram pressure of free-falling gas to the magnetic pressure of 
the field. The inner radius of the disc is thus determined by the 
magnetic field strength (assumed dipolar with a magnetic moment 
fj, = Bsr*), stellar mass and the accretion rate onto the star: rH 2 oc 
jU 2 M^~ 1/,2 rh — 1 . The inf ailing material will t hen add its angular mo- 
mentu m to the star, causing it to spin faster. Illlarionov & Sunvaevl 
( 1975) noted that if the magnetospheric radius extends beyond the 
corotation radius (r c = (GM, ) 1//3 f27 3,/2 ; where the Keplerian fre- 
quency of the disc equals the star's rotation frequency), the stellar 
magnetic field will spin faster than the gas at the inner edge of the 
disc, generally called the 'propeller' regime. A standard view ex- 



pressed in the literature is that in this case the mass in the disc is 
expelled from the system (see below) and the star is spun down. The 
exact relationship between the spin evolution of the star and the ac- 
cretion rate will depend on the physics of the disc-magnetosphere 
interaction when the inner edge of the disc is close to r c . 

Although observations of magnetic accreting stars generally 
support this general picture for magnetospheric accretion, ques- 
tions remain. A number of persistent X-ray pu lsars show spin-down 
as well as spin-up (e.g. iBildsten et alj|l997l) . or a rate of spin-up 
much lower than expected based on the accretion rate. In addition, 
many sources show hysteresis, wher e the luminosity differs be - 
tween spin-up and spin-down phases dCamero-Arranz et aifeoiOl) 



or even anti-correlates with spin-down ( Chakrabarry et al. 1997). 
Finally, in persistent sources that show both spin-up and spin-down, 
the magnitude of |z> often stays nearly constant when the torque 
changes sign, which is not naturally explained in a model where 
the torque scales as a power of rh. 

More recently, two transient X-ray pulsars were observed 
to undergo brief weak outbursts (lasting a bout 6-7 days) fol- 
lowed by a short period of quiescen ce dHeinke et alj 1201 d : 
iHartman. Galloway & Chak rabartv 201(J). The outburst recurrence 
time (on the order of a month in both sources) is too sho rt to be 
natur ally ex plained by the ionization instability model (e.g. lLasot3 
1200 lh . and IHartman. Galloway & Chakrabartvl feoich have sug- 
gested (based on the total outburst luminosity) that a significant 
amount of mass could remain stored in the disc when the accretion 
rate onto the star drops by at least 3-4 orders of magnitude. This ex- 
planation is inconsistent with the standard relation between accre- 
tion rate and inner radius, in which n n oc rh~ 2 ^ 7 , so that a change 
of 10 4 in accretion rate corresponds to an increase of ~ 10 x in rj n . 

In one of these sources, NGC 6440 X-2. |Patruno et alj feOld) 
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recently reported the detection of a strong QPO at 1 Hz. A similar 
QPO was previously detected in SAX J1808.4-3658 jPatruno et ail 
120091) . This QPO period is of the order 1CT 2 - 10" 1 times the 
viscous timescale at the corotation radius r c and hints further at 
interaction between the disc and the field. 

Similar outburst time scales, in units of the viscous time scale, 
are seen in a class of young stars called 'EXors'. These stars typi- 
cally show episodic changes in luminosity (of between 2-5 magni- 
tudes) on timescales of a few years. The variability timescale im- 
plies accretion rate variability in the inner regions of the accretion 
disc, where it interacts with the magnetic field. 

Together, these puzzling observations suggest that significant 
physics is missing in current interpretations of the phenomena seen 
in accreting magnetospheric systems. 

1.1 Accretion at a centrifugal barrier 

One of the things missing in the standard interpretation is the in- 
sight that mass prevented from accreting by the presence of a 'cen- 
trifugal barrier' does not necessarily have to leave the system at all. 
In other words, mass transferred from a binary companion onto a 
spinning magnetosphere does not have to be 'propeller e d' out . 

This was alread y noted by ISunvaev & Shakurj Jl977» and 
ISpruit & Taarnl Jl993h . Without mass loss, accretion on the mag- 
netosphere can instead be cyclic, with periods of accumulation out- 
side the corotation radius alternating with accreting phases. The 
problem of disc responding to the torque exerted at its inner edge 
by a spinning magnetosphere has a well-defined solution within the 
standard thin viscous disc formalism without mass loss. Mass loss 
in a magnetically driven wind may happen as well, but this bit of 
(still poorly known) physics is separate from the effect of a cen- 
trifugal barrier on a viscous disc. 

We have studied this pro blem in iD'Angelo & Spruitl feOld) 
and lD'Angelo & Spruitl feoill) (hereafter DS10 and DS11 respec- 
tively). There we classified the properties of time dependent so- 
lutions of the thin disc diffusion equation with a magnetic torque 
acting at its inner edge (DS10), and studied the long term evolu- 
tion of discs with such torques, including the spindown of the star 
(DS10). 

Since the magnetic torques transport angular momentum out- 
ward, the density profile of the disc is altered substantially. The disc 
remains truncated only slightly outside r c even when the accretion 
rate declines by several orders of magnitude. We found that the for- 
mation of a dead disc often results in the instability described by 
Sprui t & Taarnl dl993l) , which we suggested could be operating in 
EXors as well as SAX J1808.4-3658 and NGC 6440 X-2. In the 
instability, the disc is initially truncated outside r c and accretion 
is suppressed, causing material to build up in the inner regions of 
the disc. Eventually, the surface density in the disc becomes high 
enough to overcome the centrifugal barrier and accrete onto the 
star. Once the reservoir is emptied, the disc again moves outside r c 
and the cycle begins again. 

In DS10 we concluded that the presence of the instability de- 
pended on the mean accretion rate, rh and the details of the disc- 
field interaction.We found that the uncertainties in the physics of 
the disc-field interaction could be parameterized into two indepen- 
dent, unknown (but constrained) length scales: the radial width of 
the region of interaction, Ar, and the length scale over which the 
disc moves from an accreting to a non-accreting state Ar2. We 
made a preliminary investigation of the unstable solutions for dif- 
ferent values of our independent parameters and found the shape 
and duration of the cycles of accretion varied widely. In the present 



work we investigate the parameter space more thoroughly in order 
to characterize the instability better and link it with observations. 



1.2 Trapped discs 

The modes of accretion of a viscous disc on a rotating magneto- 
sphere can be classified in three states, all of which have a different 
appearance and effect on the star's evolution. We explicitly neglect 
the possibility of outflows, so that the disc never crosses into the 
propeller regime for any rh. 

The disc state depends on the ratio of the mean accretion rate 
through the disc (m) to the accretion rate that puts n n at r c (rh c ): 

(i) rh S> rh c : ri n <r c , star spins up 

(ii) rh ~ m c : ri n «r c , star spins either up or down, inner edge 
stays near corotation while rh varies. 

(iii) rh ~ 0: ri n >r c , star spins down, negligible accretion (dead 
disc). 

State (ii) can be further divided into cases where accretion takes 
place in a continuous way, and cases where accretion is cyclic, with 
bursts of accretion alternating with quiescent phases. Our main in- 
terest in the following paper is this second case of cyclic accretion. 
If there is a long-term accretion rate imposed at the outer edge of 
the disc, as in a mass-transferring bi nary, dead disc states can oc cur 
as part of the cycle, as suggested bv lSunvaev & Shakural Jl977l) . 

In DS 1 1 , we turned our attention to the long-term evolution of 
a viscous disc, and its effect on the spin evolution of the star. This 
added a new parameter to the problem, the moment of inertia i» 
of the star, and introduced a new characteristic timescale, Tsd, the 
spin-down timescale of the star. 

When initial conditions are chosen such that r ln is initially 
inside r c , but moves gradually outward, we found that the evolu- 
tion of the disc often gets 'trapped' with ri n near corotation, with 
a slowly decreasing accretion rate. The same happens when r- m is 
initially outside r c , but spindown of the star causes r c to move out. 
When it catches up, ri n continues to hover slightly outside coro- 
tation as the two move outward. The outward angular momentum 
transport due to the magnetic torque is accompanied by low-level 
accretion onto the star. We called this phenomenon a 'trapped disc'. 
It is the intermediate state ii) mentioned above. 

The disc does not stay trapped in all cases, however: ri n can 
also evolve well beyond r c to become a 'dead disc' - the state iii) 
above. Whether or not the disc becomes trapped depends n the de- 
tails of the disc-field interaction and on the ratio T v i sc /TsD (where 
T v i sc is the viscous accretion timescale of the disc at r c ). 



1.3 Observations and Spin Evolution of Stars with Trapped 
Discs 

Accretion in the trapped state can be continuous or cyclic. In DS1 1 
we limited ourselves mostly to parameter regimes that lead to con- 
tinuous accretion, since the short time scale of the cycles makes it 
difficult to follow the long-term evolution of the disc. In this paper 
we study the cyclic case of trapped discs in more detail, with less 
emphasis on the system's long-term evolution. We study in particu- 
lar the modification to the expected relationship between spin-down 
torque and luminosity, and the appearance of cyclic accretion. 

We explore the nature and interpretation of these cycles in 
more detail here and study how the characteristics of the instability 
(such as its period and amplitude) change with the model param- 
eters. These are discussed in section [4] When accretion is cyclic, 
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the torque between disc and magnetosphere varies over a cycle; its 
average over a cycle is what determines the net spindown or spinup 
torque on the star. In sec. 14.41 we study how this net torque differs 
from the steadily accreting case. Before that, in section [3] we dis- 
cuss the observability of the quiescent "trapped disc", and the spin 
evolution of the star as a function of accretion rate. 

In Sec. [5] we discuss our results in the context of the various 
observations mentioned above (persistent X-ray pulsars, the short 
recurrent outbursts in X-ray pulsars, and the episodic accretion 
bursts seen in both X-ray pulsars and young stars). We argue there 
that trapped discs may be related to a number of currently unex- 
plained phenomena seen in magnetospherically accreting sources. 



2 THE MODEL FOR MAGNETOSPHERIC ACCRETION 

We briefly summarize the characteristics of our model for magneto- 
spheric accretion, reviewing our description of the disc-field inter- 
action and how this alters the structure of the disc when n n > r c . 
For a compact description see DS 11; for more detail and our nu- 
merical implementation see DS10. 

As mentioned above, the interaction between an accretion disc 
and magnetic field is likely confined to the innermost regions of 
the disc, so that most of the disc is shielded from the field. The 
coupling between the disc and the field distorts the field lines 
by differential rotation, which generates a toroidal field compo- 
nent and exerts a magnetic torque on the disc. In a very short 
time, the field lines become sufficiently distorted that they in- 
flate and open, which can temporarily sever the connection be- 
tween the field and disc, before reconnection events reestablish 
connection with the star. As a result, the variations in the mag- 
netic field associated with the magnetosphere interaction will take 
place on timescales of order the rotation per i od of the star. P, |Alyl 
j 19851: lHavashi, Shibata & Matsumotd[l996l : iMiller & Stondl 1997k 
iGoodson. Winglee & Boehirn]|l997h : much shorter than the cycle 
times, which take place on a variety of viscous time scales in the 
disc. 

It is sufficient to adopt a time-averaged value for the strength 
of the generated toroidal field component, 77 = B<f,/B z and the 
width of the interaction region Ar/r- m . The size of the interac- 
tion region is unknown, so it is a free parameter in our model; we 
assume only that Ar/r < 1, and explore values in the range sug- 
gested by numerical simulations. (For a more extensive discussion 
of time variability in 77 and Ar see DS10.) 

When the inner edge of the disc r- m is outside corotation r c , a 
magnetic torque Tb acts at the inner edge. The condition that this 
torque is transmitted outward by the viscous stress requires: 



exerted by the field is 



37ri/S(n n )r f2 K (n n ) = Ta- 



il) 



This yields a boundary condition on the surface density at ri n . 
When ri n is inside r c however, the viscous torque and the surface 
density at rj n vanish: the standard case of a viscous disc accreting 
on a slowly rotating object applies. The viscous torque at n n thus 
varies rapidly over the transition width Ar around r c . This can be 
taken into account in (TJ by making Tb a function of r- m , 



T B = yv(r in )T , 



(2) 



where j/e is a suitably steep function with the property that it varies 
from to 1 over the width Ar: j/s — >• for r c —r- m 3S> Ar, j/e — > 1 
for rin — r c ^> Ar. For To we approximate the field strength at 
ri n as that of the star's dipole component, so the magnetic torque 



To = r\\x Ar/r 



(3) 



where /1 is the star's dipole moment and 77 a coefficient of order 
unity. For numerical reasons, it is important that ys be a smooth 
function; we use a 1 + tanh-profile of width Ar. For the viscosity 
we use a fixed radial dependence, v = fc„r 1//2 , corresponding to an 
a-disc of constant a and aspect ratio H/r. 

When ri n < r c the disc can accrete onto the star. The position 
of the inner edge of the d isc is then related to th e accretion rate. We 
use a standard estimate (Spru it & Taaml Fl993). which equates the 
magnetic torque with the torque needed to keep the accreting mass 
corotating with the star. It can be written as 



4fi*r? 



(4) 



(the subscript a stands for 'only in the case of accretion') and holds 
only for ri n <r c . For ri n >r c , the accretion rate drops to zero due to 
the centrifugal barrier. The relation between ri n and rh thus also 
changes steeply around r c . This is another critical element that 
needs to be taken into account. We incorporate it by writing: 



rh co (rin) = JAn(rin)m a (rin) 



(5) 



where j/ m (ri n ) describes the steep variation in the transition zone, 
and rh co is the mass flux measured in frame comoving with the 
inner edge of the disc: 



m co = rh + 27rri n E(rin)rin 



(6) 



where n n is the rate of change of the inner edge radius (the dis- 
tinction between rh and 77i co was a crucial point missing in the for- 
mulation of ST93). For the 'connecting function' y m , with limiting 
values 1 inside corotation and outside, we again use a smooth 
tanh- function. The width of this transition must be of the same 
order as the transition width Ar introduced above, but the physics 
determining it is rather different, so we keep it as a separate pa- 
rameter Ar2, and explore it independently.By introducing Ar and 
Ar2 as the transition length scales from one boundary condition 
to other, we reduce the uncertainties of the transition region to two 
free but constrained parameters, which we vary systematically to 
explore the range of unstable solutions. 

Eq. © and ^ define an equation for ri n (t). Together with the 
standard thin disc diffusion equation, we have a set of equations, 
second order in time and space, for £ and ri n . The equations are 
stiff inside corotation since one of the time derivatives (for n n ) 
disappears. A numerical method appropriate for this circumstance 
has to be used. 

The steady-state case, dt = 0, is illustrative. The solution of 
the thin disc diffusion equation can be written as 



37TZ/E = 



T B 



1/2 

— ) +rh 



1- — 



1/2' 



(7) 



n(n»)r£, 

where Tb is given by (O and 771 by l[4][6j. 

For rin well outside corotation (ri n - r c S> Ar), 771 vanishes 
and we have a dead disc. The surface density is then determined 
by the first term on the RHS. The steady outward flux of angular 
momentum in this case has to be taken up by a sink at some larger 
distance, otherwise the disc could not be stationary as assumed. 
This sink can be the orbital angular momentum of a companion 
star, or the disc can be approximated as infinite. The latter is a good 
approximation for changes in the inner regions of the disc, if we 
consider short time scales compared with the viscous evolution of 
the outer disc. 
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2.1 Characteristic numbers 

We choose r c (the star's initial co-rotation radius) for our length 
scale, and define t c = vjr 2 , the viscous accretion timescale at 
r c as our characteristic timescale. We adopt vq = a(H/r) 2 = 
10 -3 for the dimensionless viscosity parameter, and r\ = 0.1 as the 
relative size of the component. We also define a characteristic 
accretion rate m c : 

*■ = 

This is the accretion rate in {4]l that would put the magnetospheric 
radius at r c . 

For a typical ms-pulsar: 

-5/3 / „ \ 2 



m c = KT^Moyr" 1 



1.4M, 



V10 8 G J 



(9) 



/_B. 



V10 e cm/ ViOOms 
while for the case of a T Tauri star disc, 



-7/3 



m c = 1.4 x lO~ 9 M yr _1 



M. 
0.5M 



-5/3 



-7/3 



(10) 



V1500Gy V 2i W V 7d 

To characterize the spin change of the star, we adopt a characteristic 
torque based on the definition for m c : 

j c =rh c (GM,r c ) 1/2 , (11) 

which is the rate at which angular momentum is added by mass 
accreting at m c . This definition is somewhat arbitrary, but has the 
advantage of being independent of Ar/r- m and Ar2/nn. 



3 SPIN EVOLUTION AND PHYSICAL PROPERTIES OF 
A TRAPPED DISC 

3.1 Observability of trapped discs 

The density structure of a trapped disc deviates from the standard 
accreting disc. In a dead disc (the limiting case of a trapped disc), 
the accretion rate is zero and the surface density profile is deter- 
mined by the rate of angular momentum transport from the star to 
the disc. 

Provided that the thin disc approximation remains a good de- 
scription for the disc, meaning that the energy generated by viscous 
turbulence is dissipated locally the temp erature of the dead disc can 
be estimated dSunvaev & Shakurall 19771) from the viscous dissipa- 
tion. 

The local energy dissipation rate in a viscous disc is: 

Q + = uT,{rQ' K ) 2 = ^Efi|. (12) 

Near the inner edge of the disc where most of the energy is dissi- 
pated, the disc can be approximated as steady, so Q applies. For a 
dead disc (rh = and n n well outside corotation), Tb = To, and 
the surface density as a function of r is given by 



= — ( To 
3rr V Mr 2 



If this energy is radiated locally as a blackbody, then 

Q + = Q- = 2aT* 



(13) 



(14) 



where T s is the surface temperature of the disc as a function of 
radius (the factor 2 arising because the energy is lost from the two 
surfaces of the disc), and <tb the Stefan-Boltzmann constant. With 
< f3TT3l this yields 



<tbT s 4 



3 2 



Ar 



GM 



1/2 



(15) 



which varies with r as r 7 ' 2 , showing that most of the energy is 
radiated away near the inner edge by the mass piled up at the 'cen- 



,1/2 



as we assume in the calculations to 



trifugal barrier'. If v r. 
follow, J 1 5 b implies T B ~ r _1 . 

With the nominal object parameters used in 12. II a dead disc 
with its inner edge near corotation will have a maximum temper- 
ature, at ri n , of T s ~ 3007^ in the case of a T Tauri star and 
T s ~ 20000 K for a ms X-ray pulsar. This estimate, however, only 
includes the internal energy dissipation in the disc. Both in proto- 
stars and X-ray binaries, reprocessing of radiation from the central 
object by the disc usually dominates over internal dissipation, espe- 
cially at larger distances from the centre. Nonetheless, it might be 
possible to identify a disc as belonging to the dead class with more 
detailed information on the spectral energy distribution. 



3.2 Accretion rate and angular momentum exchange with 
the star 

The interaction between the magnetic field and the disc outside r c 
removes angular momentum from the star, while the gas accreting 
onto the star spins it up. The net torque on the star is thus a function 
of both the average accretion rate through the disc and the location 
of ri n . It is given by: 



i = m(GM trm ) 1/2 -^^!/ S , 



(16) 



where the first term gives the spin-up from accretion and the second 
term is the spin-down from the disc-field interaction. (The spin- 
down torque vanishes for r- m < r c , which we impose with the j/e 
smoothing function.) In sec. [2] we defined the relationship between 
m and ri n , so that J will just be a function of m, Ar/r{ u and 
Ar2/Vi n . In the rest of the section we illustrate the implications of 
this relationship for the spin-evolution of the star. 

Fig. Q] shows the relationship between ri n and m as Ar2 /ri n 
changes, to demonstrate how the disc can become trapped. The 
three curves show r- m as a function of rh for Ar2/r- m = 0.05 (solid 
blue), 0.1 (dashed green) and 0.5 (triple-dot dashed red). At high 
accretion rates (ri n < r c ), all three curves scale as ri n oc m~ 1//5 . 
However as rh decreases the solutions quickly diverge. For small 
Arz/r-m (Ar2/Vin =[0.05, 0.1]), there is a knee in the solution 
around r- m = r c , and r- m increases much more slowly as m de- 
creases, so that ri n remains close to r c even when rh has decreased 
by several orders of magnitude. For Ar2/r- ln — 0.5, the knee 
straightens out since the transition from accreting to non-accreting 
solutions is much more gradual. 

Fig. [2] shows the dependence of the net torque J on the star 
on accretion rate rh, for different values of [Ar, Ar2]. At high ac- 
cretion rates, the disc is strongly accreting, star spins up with the 
torque j oc m 9 ^ 10 . At lower accretion rates, the interaction region 
moves outside r c , causing the net angular momentum exchange to 
change sign, and spinning the star down. As the accretion rate de- 
creases further and n n moves increasingly far outside r c into the 
effectively 'dead' state, the spindown torque decreases because of 
the decline of the star's field strength with distance. 
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Figure 1. Change in r- ln as a function of m for different Ar2/fi n . Solid 
curve: Ar2/r; n = 0.05, Dashed curve: Ar2/Vj n = 0.1, Triple-dot-dashed 
curve: Ar2/r; n = 0.5 
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Figure 2. Change in net torque on the star as a function of rh for different 
values of Ar/n n and Ar2/n n . From top to bottom, the individual panels 
correspond to Ar2/r; n = [0.05,0.1,0.5]. The different patterned lines 
correspond to different values of Ar/r- In : 0.05 (solid), 0.1 (dashed), 0.5 
(dot-dashed), and 1.0 (triple-dot dashed). 



The detailed shape of the curves in figs. UI2l reflects the (tanh- 
) shape of the transition functions ys and y m we have used. More 
generally, however, when J < 0, the amount of spin down is di- 
rectly proportional to Ar/n n , and Ar2/n n determines how far 
away from r c the disc can move as rh decreases. 



4 CYCLIC ACCRETION 

The possibility of identifying cycles like those seen in DS10 in ac- 
tual observations is a interesting prospect. Cycles would provide 
a direct observational connection with the trapped disc state phe- 
nomenon that we have identified here as a likely consequence of 
the disc-magnetosphere interaction. The shape of the cycles and 
circumstances of their occurrence would provide important clues 
about the details of the interaction region, which we have simply 



parametrized with the two transition widths Ar/r m (far the torque) 
and Ar2 /r- m (for the mass flux). 

In this section we expand the analysis in DS10 to investigate 
the properties of the instability more quantitatively. In Sec. 14. II we 
investigate the presence of the instability as a function of the sys- 
tem parameters. In sec. 14.21 we show how the period and amplitude 
vary with m, Ar/n n and Ar2/n n . In sec. 14.41 we investigate how 
the presence of cycles will change the spin evolution of the star in 
comparison to a steadily accreting disc. Finally, in sec. 14.51 we ex- 
amine the appearance of cycles for non-steady accretion in systems 
in which the star's spin is evolving. 

4.1 Parameter map of the instability 

In DS10 we showed that the disc instability depends on Ar/n n , 
Ar2/n n and rh, and ran simulations in the parameter spaces 
[Ar, rh] and [Ar2, rh] to determine when the instability occurred. 
Additional simulations of trapped discs at very low accretion rates 
show the instability is present for a larger parameter space than we 
investigated in DS10, so here we have repeated our analysis on a 
larger parameter space and better-sampled grid. 

As in DS10, we keep the rotation rate of the star fixed, since 
the time scale of the cycles is much shorter than the time scale of 
spin changes of the star. The corotation radius is then a fixed dis- 
tance, and defines a unit of mass flux, m c given in sec. 12.11 As 
unit of length we can use r c , and as unit of time the viscous time 
scale r^/u(r c ). Apart from a parameter specifying the radial de- 
pendence of v, which we keep fixed throughout, the problem is then 
defined by the three dimensionless parameters Ar/rin, Ar2/n n , 
and m/m c . 

It turns out (discussed below) that this 3-dimensional param- 
eter space contains 2 nearly separate regions of instability, which 
can be characterized using a few 2-dimensional slices instead of 
having to scan the entire space. 

In comparison to DS10, the present results survey a larger pa- 
rameter range at higher resolution for all three variables. The in- 
dividual simulations at higher spatial and time resolution. In total 
we ran 1545 simulations: 855 in a 2-dimensional slice [Ar, rh], and 
690 in a [Ar2, rh] slice. Our goal was to determine the regions in 
parameter space where cyclic behaviour occurs, as well as how the 
characteristics of the cycles (such as period, amplitude and appear- 
ance) change as a function of the parameters. 

Figures[3]and|4]show the unstable regions across [Ar, rh] and 
[Ar2,rh]. In fig. [3] we kept Ar2/n n fixed at 0.014 and surveyed 
the range Ar/r in = [0.02,0.28] and rh/rh c = [10" 6 , 6], Fig- 
ure [4] surveys the same range in rh and the range Ar2 /ri n 
0.00 1 , 0.14], keeping Ar/ri n fixed at 0.05. The figures (especially 
figure [4]l show that there are two distinct regions in which the in- 
stability is active, which we call RI and RII. 

The instability has different properties in the two regions: the 
period, shape and amplitude of the outburst are all qualitatively dif- 
ferent, as we demonstrate further below. 

The first instability region (RI) appears in a small range 
of Ar2/ri n = [0.002,0.03], but over a considerable range in 
Ar/ri n (up to 0.25) and five orders of magnitude in accretion rate: 
rh/rh c ~ [10 -6 , 10" 1 ]. This instability region was the focus of our 
study in DS10, where we discussed in detail the appearance of the 
instability. (The phenomena summarized below can be seen in figs. 
4,5 and 6 of that paper.) The instability is characterized by large am- 
plitude outbursts followed by long periods of quiescence in which 
the accretion rate onto the star drops to zero. The duty cycle for 
the instability depends on all three parameters (Ar/ri n ,Ar2/ri n , 
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from the output of our [Ar2 , rh] set of simulations. The simula- 
tions all have m = 1.05 rh c and Ar/n n = 0.05, with increasing 
Ar 2 /r in from bottom to top: Ar 2 /r in = [0.028, 0.05, 0.1, 0.15]. 
The instability has a different character from the RI instability. The 
quiescent phase is absent and there is always some accretion onto 
the star, even during the low phase of the cycle. The shape of the in- 
stability is also significantly different from the RI instability: there 
is no additional higher frequency oscillation, and the outburst pro- 
file is smoother. In particular, the initial spike of accretion seen in 
RI is completely absent. The profile is nearly sinusoidal for small 
values of Ar^/nn (although the duty-cycle is always less than 0.5). 
For larger values of Ars/nn, the outburst is characterized by a 
rapid rise, followed by a decaying plateau and then rapid decline to 
the low phase of the cycle. 



Figure 3. Instability map in [rh, Ar/rjJ, keeping Ar-2/Vin = 0.014 
fixed. The shaded regions denote unstable simulations. The instability oc- 
curs in two nearly disconnected regions of the parameter space. In this slice, 
one (RI) is present over a large range in Ar / r; n and rh, the other in a small 
region in Ar/ri n at rh ~ rh c (RII). 




0.02 0.04 0.06 0.08 0.10 0.12 0.14 

Ar 2 /r in 



Figure 4. Instability map in [rh, Ara/rjn], keeping Ar/r- ln fixed at 0.05, 
with shaded areas denoting unstable cases. The two regions of the instability 
are more clearly separated than in the slice of fig. [3 RI only occurs for very 
small values of Ar^/r^ while RII occurs around rh ~ rh c , over a large 
range in Ar 2 /ri n , extending beyond the region of the graph. 

and rh), decreasing significantly as the mean accretion rate drops 
(and mass takes longer to accumulate and power another cycle). 
The outburst profile generally takes the shape of a relaxation os- 
cillator, with an initial accretion peak followed by a tail of much 
lower amplitude accretion. Additionally, during the long phase of 
the outburst suboscillations sometimes appear in the accretion pro- 
file. These suboscillations have a much higher frequency than the 
overall burst, and appear to be the second instability (RII) super- 
imposed on the outburst while the mean acccretion rate onto the 
star is temporarily higher. As the mean accretion rate is decreased, 
the outburst becomes shorter and shorter, until finally it simply ap- 
pears as a single spike of accretion followed by a long period of 
quiescence. 

The second instability region (RII) occurs only around 
rh/rh c ~ 1, in a relatively small range of Ar/n n = [0.01, 0.07], 
but over a large range in Ar2/n n . (We have truncated the figure 
at Ar2 /Yin to make RI more prominent, but additional simulations 
show that the instability continues to larger values of Ar2/ri n , up 
to at least 0.4). Fig [5] shows four sample simulations for RII taken 



4.2 Period and Amplitude of Instability 

We can use our high resolution survey of the instability parameter 
space to study quantitatively how the properties of the instability 
change as a function of the different parameters. This allows us to 
see the differences between the RI and RII instability regions more 
clearly, and better understand how the instability operates. It also 
allows us to make more general predictions for the appearance of 
the instability that can be tested against observations. 

To measure the period of the instability we autocorrelate the 
output of each simulation, and take the first peak in the autocor- 
relation. All our simulations take as initial conditions the stable 
solution given by (|7) that then becomes unstable. However, the in- 
stability generally takes some time to reach a stable period and am- 
plitude, which varies significantly between simulations. This intro- 
duces some error in the estimate of the period (particularly in RII, 
where the instability takes longer to emerge), which is reflected in 
the figures below. 

Figs. [6] and [7J show the instability period [top] and maximum 
amplitude [bottom] as a function of accretion rate for each unstable 
simulation in figs.[3]and|4]respectively. Each individual curve rep- 
resents a different value for Ar/ri n or Ar<2/r- m , with the colours 
evolving from purple (smallest parameter value) to red (largest). In 
the top panels, the gaps in the curves denote stable cases. In the 
bottom panels, the dashed y = x line shows the mean accretion 
rate for each simulation. 

Some general properties of the period and amplitude of the 
outburst are valid for both for RI (low rh; cf. sec. 14. It and RII (high 
rh). The period of the instability shows an approximately power- 
law dependence on accretion rate, with an index independent of 
Ar/rin and Art/ Yin. It also varies strongly with both Ar/n n and 
Ar2/n n . The period scales directly with Ar/nn. The maximum 
amplitude of the cycle also scales roughly linearly with Ar/nn 
andAr2/ri n . It shows the opposite trend for Ar 2 /ri n , decreasing 
as Ar 2 /ri n increases. On the other hand, the amplitude of the out- 
burst is nearly independent of the mean accretion rate onto the star 
for a given Ar/nn or Arz/nn, as is seen in the bottom panel of 
figs. [6] and [JJ Rather than producing a weaker outburst, at low rh 
the duration of the low phase of the outburst increases as mass ac- 
cumulates in the disc. 

This result demonstrates that the details of the outburst, apart 
from the period of the cycle, are largely independent of rh, and 
depend almost only on Ar/n n and Ar^/nn- For large Ar/n n a 
larger reservoir of mass is available for outburst, so that the dura- 
tion of each burst increases. Smaller Ar2 /n n means a more abrupt 
transition around r c from accreting to non-accreting states, so that 
for the same rh, a simulation with small Ar2 /n n can build up more 
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Figure 5. The accretion instability in RII, the second instability region for Ar/n n = 0.05 and fa = 1.05 fa c , with Ar2/ri n = [0.028, 0.05, 0.1, 0.15]. In 
each figure, the solid line shows the evolution of r; n over a cycle, while the dot-dashed line shows the accretion rate scaled to units of fa c . The dashed line 
shows the inner radius corresponding to steady-state accretion. The instability only manifests around fa ~ fa c , but occurs over a large range in Ar'2/r; n , the 
transition length between accreting and non-accreting states. Unlike the RI instability, the accretion rate never drops to zero during the quiescent phase of the 
instability. 



mass in quiescence than a simulation with larger Ar2/ri n , leading 
to a longer outburst phase and longer period. 

Although the trends described above apply to both RI and RII, 
figs. [6] and |7j clearly demonstrate significant differences between 
the two regions. In RI the accretion amplitude scales inversely with 
both Ar /n n and Ar2 / r- m , so that a large Ar /r- m (or Ar2 /n n ) cor- 
responds to an outburst with an amplitude up to 100 x smaller than 
for the same fa and smaller Ar/ri n (Ar^/fin). The transition be- 
tween accreting and non-accreting states becomes more gradual as 
Ar/nn (Ar2/n n ) increase, so it makes sense that the initial peak 
of the outburst will be lower for larger Ar/ri n (Ar2/ri n ). Interest- 
ingly, the cycle amplitude scales inversely with Ar/rj n , while the 
period scales directly. The opposite behaviour is true of Ar2/ri n . 
Instabilities with large Ar/r- m thus have weaker outbursts and a 
longer duty cycle (compared to an instability with smaller Ar/rin), 
while those with a large Ar^/r^ instead have weaker outbursts 
and a shorter period, and less variation between outburst and qui- 
escence. 

The RII instability is presently not as well-sampled as RI, but 
some general characteristics are evident. Most significantly, the am- 
plitude is much larger (~ 40%) than the largest amplitude RI in- 
stability, and the period of oscillation is up to ~ 60% shorter. As 
well, our present results suggest that the outburst amplitude for RII 
is roughly independent of Ar/ri n , Ar2/ri n , and m, staying fixed 
at about 10 rh/m c . Instead of changing in amplitude, as (m) de- 
creases (or Ar2/n n increases), the period increases as well and the 
source spends more time in the low phase of the cycle (as is seen in 

fig.EJ. 

In summary, the instability is very different in RI and RII. In 
RI, the period of the instability is 1-1000 times longer than the 
viscous accretion timescale in the inner parts of the disc, and is 
a strong function of Ar/rj n , Ar2/ri n and fa. The amplitude and 



duty cycle of the outburst are strong functions of the detailed disc- 
field interaction around r c , but appear nearly independent of m. 
The instability occurs over a wide range in fa and Ar/ri n , but is 
confined to small values of Ar2/n n , the parameter describing the 
transition from accretion to 'centrifugal barrier'. RII is confined to 
a small range in fn, but extends to much larger values of Ar2/n n . 
The cycle period is typically shorter than the viscous timescales 
in the inner disc (P cyc ~ 10 -2 — 1 t c ), and the amplitude of the 
outburst is higher than RII. The amplitude of the outburst appears 
to be independent of fa, Ar and Ar2, with a typical value of ~ 
10 fa c . 

4.3 Interpretation of the instability regions 

The instability in region I is easie st to understand; it co r respo nds to 
the cyclic accretion predicted b y ISunvaev & Shakural Jl977l) . and 
studied in Sprui t & Taariil (l993). Mass piling up at the centrifugal 
barrier eventually becomes large enough to 'open the gate', and the 
resulting accretion rate is large enough to push the inner edge of 
the disc inside corotation for a while, until the mass reservoir has 
been drained. Unlike a fixed reservoir, the mass in it is not a fixed 
number, since it depends on the extent of the disc participating in 
the cycle. This explains the large range in cycle period over which 
it operates dSprait & Taamll 19931) . 

This view of the RI cycle implicitly assumes a sharp transition 
from accretion to pile-up outside r c . If the transition is not sharp 
and the centrifugal barrier is 'leaky', the cyclic behavior can be 
avoided. Mass still piles up outside r c , but accretion can be matched 
by mass leaking through the barrier. This e xplains why the cycle s 
are limited to low values of Ar2 / r- m . As in lSpruit & Taaml dl993h . 
the cycles can start from a steady accreting state in the form a linear 
instability, which sets in only if the transition is steep enough. 
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Figure 6. Period [top] and amplitude [bottom] of the instability as a func- 
tion of rh I rh c for the unstable simulations shown in fig. [3] The accretion 
rate spans the range rh = [10" 6 ,6] and Ar/r in = [0.02,0.28]. The 
individual curves show differen values of Ar/rj n , increasing from black 
(smallest) to red (largest). In the top panel, the gaps in the curves indicate 
regions of stable accretion. In the bottom panel, the mean accretion rate is 
indicated by the dashed y = x line. 



Region II is restricted to accretion rates around the transition 
to steady accretion , and small values of Ar/Vi n . Its restricted range 
of oscillation periods - a fraction the viscous time scale at r c - 
shows that it operates on a region of finite extent, of the order of 
the width of the transitions Ar/n n , Ar2/ri n . In the transition re- 
gion, a small increase in accretion rate will cause n n to move in- 
wards, which will substantially decrease the torque at n n (because 
Ar/rin is small) and moves n n even further inward. As n n moves 
inward the accretion rate through n n increases and an outburst is 
triggered. Since even in the low phase of the cycle material is ac- 
creted onto the star, the reservoir of matter available for the outburst 
is small and only the innermost regions of the disc are involved in 
the instability (of order Ar2/ri n ). The timescale for the instability 
is thus much shorter than the RI instability, and the period scales 
with Ar-z/nn. 

Since the instability occurs over a small region around r c , it is 
likely sensitive to the details of the physics of the transition region. 
As long as we may assume the transition to be monotonic, however, 
it's shape is rather constrained by its values far from corotation 
(0 and 1 respectively), and its slope at corotation, (controlled by 
the width parameter, Ar2). While the results above show that the 
existence and properties of the cycle depend strongly on Ar2, one 
might wonder to what extent the results also depend on remaining 
degrees of freedom in the choice of transition function. To test this, 



Figure 7. Period [top] and amplitude [bottom] of the instability as a func- 
tion of rh/rh c for the simulations shown in fig. [4] The plot is the same 
as fig. [6] replacing Ar/r; n with Ar2/r; n . Each curve shows a different 
Ar2/V; n spanning the range Ar2/r; n = [0.004, 0.15]. 



we re-ran a series of simulations using the error function rather than 
the tanh — function. Since the width of the function is a matter of 
definition in both cases, we use slope of the function at corotation as 
a well-defined measure of width for the case-by case comparison. 
The range of unstable values in the parameter space, the period, 
shape, and amplitude of the cycle were then found to be very similar 
as when using a tanh profile. 



4.4 The effect of cycles on angular momentum exchange 

The accretion instability can have a significant effect on the spin 
change with the star. When accretion is stable, the rate of angular 
momentum exchange between the star and the disc is given by H6\ , 
and illustrated in fig. [2] When accretion proceeds via cycles, the 
amplitude of the outburst is generally much higher than the mean 
rh, and the disc spends a significant amount of time in the low phase 
of the cycle, in which little or no accretion takes place (sec. 14.21 . 
This will affect the net angular momentum exchange between the 
star and the disc: since the low phase of the instability almost al- 
ways lasts longer than the outburst phase, the star will spin down 
faster than if accretion proceeded steadily. 

Using the results of the simulations presented in sec. 14. 11 and 
14.21 we can compare the torque exerted on the star for stable and 
unstable accretion: 



AJ 



(<-^sim) Jet. 
J cons 



(17) 
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Figure 8. Effect of cycles on the average torque, as a function of accre- 
tion rate. Figure shows the difference between the torque in the unstable 
initial state before the cycle develops, and the value measured in the fully 
developed limit cycle. The different curves represent different Ar/r- ln , in- 
creasing from purple (Ar = 0.05) to red (Ar = 0.22). Ar2 = 0.014 is 
held constant. Cases that do not become cyclic lie on the red line along the 
x-axis. 
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Figure 9. Difference between the torque in the unstable initial state be- 
fore the cycle develops, and the value measured in the fully developed limit 
cycle. The different curves represent different Ar2/r ln , increasing from 
purple (Ar2 = 0.004) to red (Ar = 0.028). Ar = 0.05 is held constant. 
Solutions that lie along the x-axis do not show cyclic variations. 



Here J CO ns is the torque exerted on the star in the initial state (cal- 
culated using|16t, and { J s i m ) is the average torque exerted over one 
full cycle once a stable limit cycle is reached. 

Figs. [8] and [9] show AJ/J CO ns as a function of accretion 
rate for the simulations in figs. [3] (varying Ar/n n ) and [4] (vary- 
ing Ar2/n n ), respectively. The individual curves show different 
Ar/rin (or Ar2/ri n ), with the colours from purple to red indicat- 
ing increasing values of the parameter. Lines along the x-axis are 
stable solutions. 

The results are easily separated into the RI (low rh) and RII 
(high rh) instability regions. In RI, both J CO ns and J a i m are neg- 
ative (cf. [2} over the entire unstable region, and the star spins 
down. As expected, the presence of the instability increases the ef- 
ficiency of the spin-down torque considerably: by up to 50% for 
[Ar/Vin, Ar2/)"in] = [0.05, 0.015]. For a given value of Ar/ri n 
or Ar2 /ri n , A J/ J cons is largest at the highest rh and decreases for 
smaller rh. When rh stays fixed, A J/ J cons is largest for the small- 
est values of Ar/r^ and Ar2/ri n (where the outburst amplitude 




Time [T >is J 



Figure 10. Unstable disc solution for [Ar, Ar2] = [0.2, 0.02]r in . At 
t = 0, ri n /r c = 1.3, so that rh ~ 0. Initially r- m remains constant, 
before decreasing sharply as r c increases. This prompts the onset of the ac- 
cretion instability, which continues even while the star spins down. At later 
times, the amplitude and frequency of the instability decrease and the disc 
reverts back to accreting steadily. Top: Evolution in r- m (black curve) and r c 
(red curve), before, during, and after the instability appears. Bottom: Cor- 
responding accretion rate onto the star. The accretion rate is given in terms 
of m c corresponding to r c at (t = 0) . 



is largest, cf. sec. 14. 2\ . This behaviour is readily understandable: 
the added angular momentum from accretion is most relevant at 
high rh, and the rate at which this material is accreted (i.e. the duty 
cycle of the outburst) will determine how much the star will spin 
up or down. The larger the amplitude of the outburst, the faster the 
reservoir of mass will be accreted, and so the shorter the duty cycle. 

In RII the situation is quite different. In the parameter range 
for Ar/ri n and Ar2/ri n where the solution is unstable, (J s i m ) and 
Jcons > when m — rh c (fig. |2). The presence of the instability 
decreases the rate of spin-up by a maximum by up to about 10%. 



4.5 Transient instability cycles 

In DS 1 1 we studied how discs in our model evolved as the spin 
of the star changed. To do this we introduced the star's moment of 
inertia, 7» as an additional parameter of the problem, which intro- 
duced a new characteristic timescale for the problem, the spindown 
timescale Tsd'- 



T SD = P. /P. 



IM*rf n 
r\p? Ar 



(18) 



We found that the ratio of T v i sc , the viscous accretion timescale 
at r c , and Tsd plays an important role in determining whether a 
disc will become trapped with n n close to r c when rh, the mean 
accretion rate through the disc drops to zero. 

The spin change of the star can also determine whether or not 
the accretion instability occurs. For a given set of [Ar, Ar2] the 
instability frequently appears over a narrow range in rh/rh c , so 
that changing r c (and hence rh c ) will cause the instability to ap- 
pear or disappear. This effect is difficult to observe in our simula- 
tions, since the timescale of the instability (~ 10~ 2 — 10 3 T v i S c) 
is in general much shorter than the spin-down timescale of the star 
(T SD ~ 10 3 - 10 17 T visc ; DS1 1), so once the instability occurs the 
simulation typically proceeds on the instability timescale. Despite 
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this, we do find some solutions that show transient oscillations (typ- 
ically those with the long periods, P ~ 10 3 T v i sc ) for a constant m, 
which appear and disappear as r c changes. 

This transient instability is illustrated in figs. [10] and QT] For 
these simulations we have adopted the same scaling parameters 
as we used in DS11. We take r», the star's radius, as our char- 
acteristic length scale, and define t* = vfr%, the nominal vis- 
cous timescale at the star's radius, as the characteristic timescale. 
We set Tsd = 380 T v i sc (where T v i sc is the viscous timescale for 
r c (t = 0)), which is appropriate for a protostellar system. We set 
r c {t = 0) = 1.8r». 

In fig. [10] we show the evolution of an initially dead disc 
(m = 0) with [Ar/ttn., Ar 2 /V in ] = [0.2,0.02]. The disc is ini- 
tially in the steady-state solution given by Q with r- m = 1.3r c . 
The top panel shows the evolution of ri n (black) and r c (red) in 
time, while the bottom panel shows the accretion rate onto the star, 
scaled to m c at t = 0. At early times [not shown] there is no ac- 
cretion onto the star and n n remains fixed while r c moves steadily 
outwards. Once r c moves close enough to r ln (at ~ 10 4 £*) to be- 
gin accreting matter, gradients in the disc's surface density cause 
rin to move quickly closer to r c . At nn/V c ~ 1.1, however, the 
disc suddenly becomes unstable and accretion proceeds via a series 
of large amplitude bursts, which increase the maximum accretion 
rate onto the star enormously. r c continues to move outwards as the 
source oscillates, which causes m c (r c ) to decrease as well. Even- 
tually m/m c moves out of the unstable range, and the disc settles 
back into a steadily accreting solution with a gradually decreasing 
accretion rate. This result emphasizes the fact the instability can de- 
pend sensitively on m/m c , so that as r c evolves, cycles can appear 
and then disappear. 

For constant Ar/r- m and Ara/Vin, the appearance of cycles 
depends on rh and m c . As we showed in DS11, when the mean 
accretion rate through the disc drops so low that the disc becomes 
trapped, the accretion rate at ri n depends strongly on T v i sc /Tsu- 
This ratio also plays an important role in determining the presence 
and duration of unstable accretion phases. 

We illustrate this in fig.[UJ which shows a disc evolving from 
an accreting disc to a trapped disc for different T v i sc /TsD, adopting 
parameters appropriate for a neutron star. We take the field strength 
10 12 G, initial spin period 5s, and initial inner edge radius 0.95r c . 
The disc is thus initially in an accreting state. The outer boundary 
is located at 100 n n . The initial spindown timescale of the star is 
Tsd = 10 5 years. We adopt the same system parameters as for fig. 
[TOl with [Ar/r in , Ar 2 /r in ] = [0.2, 0.02]. At t = we set m = 
at the outer boundary of the disc, so that the accretion rate through 
the disc gradually drops. 

The figure shows the evolution of r c (dashed curve) and 
rin (solid curve) as a function of time for increasing r v i sc , with 
r visc /r S D = 2.5 x [10~ 9 , 10" 7 , 10~ 5 , 10" 4 ] from top to bottom. 
In all four simulations, the disc initially move outwards as rh de- 
creases to about 1.05 r c , where the accretion rate triggers the onset 
of the instability. 

The appearance of the disc changes dramatically with chang- 
ing viscosity. Discs with a high viscosity are able to diffuse out- 
wards quickly as rh declines, so that r- m in the top panel (with the 
highest viscosity) moves far from r c before the star can begin to 
spin down, and the disc never becomes trapped close to r c . The 
instability (which only occurs between m/m c — [10 -3 , 0.1]) will 
only appear briefly before ri n moves too far from r c and the accre- 
tion rate drops. In the middle two panels the instability also appears 
briefly, but the star does begin to spin down quickly enough that r c 
begins to move outwards at the same rate as n n and the disc re- 
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Figure 12. The 2008 outburst of EX Lupi. From January to August, the 
source was up to 4 magnitudes larger than its quiescent value. Figure cour- 
tesy A. Juhasz. 



mains trapped. In the bottom panel (with the lowest viscosity), the 
disc diffuses outwards slowly, and the star spins down more rapidly. 
Once the instability is triggered, n n remains close enough to r c that 
the instability persists as the star continues to spin down. 



5 RELEVANCE FOR ASTROPHYSICAL SOURCES 

5.1 Trapped discs and the disc instability in EXors 

We have previously suggested that the disc instability discussed in 
this paper is operating in a class of T Tauri stars known as 'EXors'. 
EXors, like their prototype, EX Lupi are characterized by repeated 
large outbursts: changes by up to four magnitudes in luminosity 
lastin g several months, wi th a characteristic total period of several 
years (Herbi3 l2007l l2008). Observations of these sources in both 
outburst and quiescence can produce valuable constraints on both 
the disc instability and the trapped disc itself. 

Figure[T2lshows the my brightness of EX Lupi over 500 days, 
taken from the AAVSO and ASAS archives (from Attila Juhasz, 
private communication). In January 2008, EX Lupi underwent a 
large outbu rst, increasing by 3 magnitudes, from my — 11 to 
m v = 7.9 dAspin et alJ l ioiOl) in about 20 days. In July 2008 EX 
Lupi was at my = 12.4 i Sipos et al]|2009l) . so the total change in 
luminosity was more than 4 magnitudes. EX Lupi remained in out- 
burst for eight months, decaying from maximum to about mv = 10 
before abruptly dropping to mv = 14 over about 40 days. The data 
show evidence of an additional higher-frequency modulation with 
a period of about 30 days and amplitude of 1-2 magnitudes. 

The outburst profile looks similar to some of the pro- 
files seen in our simulations, with a rapid rise and a grad- 
ual decreasing plateau phase, followed by a rapid decay as the 
source returns to quiescence. The higher frequency modulation 
is also suggestive of the higher frequency oscillations seen in 
the RI instability discussed above and in DS10, although its 
short timescale ~ 30 days, or ~ 4 — 10 x longer than P* 
dSipos et alj|2009l) might also indica te direct modulation between 

Goodson, Winglee & Boehml 



the disc and magnetic field (e.g. 



ll997llRomanova. Kulkarni & Lovelacell2008l) . 

Spectral observations taken both during the outb urst and qui- 
escence allow additional constraints on the disc model. IXspin et all 
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Figure 11. Evolution from an accreting to a non-accreting disc, for increasing (top to bottom) ratios of T v j ac /Tsu for stable discs ( Ar/ r ln = 0.1, Ar2 / r- m 

■ From top to bottom, the ratio T v ; sc /TgD is 2.5 X [10 -9 , 10~ 7 , 10 — 5 , 10 — 4 ]. The instability appears in all four panels, for varying lengths of time. 
As the ratio between the two timescales decreases, the disc is not able to move outwards as quickly before the star begins to spin down, so that rj n will always 
remain close to r c and the instability persists for longer. 



d20ld) analysed a series of spectra of EX Lupi from the 2008 out- 
burst to constrain the parameters of the disc. They calculated an 
accretion rate of 2 x lO~ 7 M0yr _1 in outburst, and modelled the 
2.2935pm CO overtone emission to suggest an inner gap of about 
10 r\ in the disc. They also used previously publishe d results to 
calcul ate 6 x 10~ 9 Mqyt^ 1 for the quiescent emission. ISipos et"ai1 
(2009) analysed a series of spectra for EX Lupi in quiescence (with 
my ~ 12), and did a detailed fit for the spectrum. They estimated 
an accretion rate of 4 x lO~ lo M0yr -1 , and fit the spectrum with a 
dust-dominated disc truncated at about 0.2 AU. These observations 
thus suggest that the accretion rate increased by ~ 30 — 500 x dur- 
ing the outburst, and that close to the star, gas is either absent or 
optically thin. 

To see whether these observations are compatible with our 
model of a trapp ed disc, we adopt the best-fit parameters of 
ISipos et alj fe0091) . with M» = O.6M , r» = 1.6-R Q , and P» = 
6.3d (based on usini = AAkm/s and an inclination of 20 de- 
grees). The magnetic field is unknown, so we assume a strong dipo- 
lar field with B, = 1500G, which is comp arable to field str engths 
observed in T Tauri stars such as BP Tau ( Donati et al . 2008). The 
rest of the parameters we adopt are the same as for our represen- 
tative model (sec. 12. lb . This puts the co-rotation radius at 7.6 r*, 
and sets ri n ~ r c for the quiescent disc. The viscous timescale at 
r c is thus ~ 30 yrs. If we assume that the period of the instability 
is ~ 60 yrs (the time between the 2008 outburst and the similar 
amplitude outburst in 1955), or ~ T v i sc . The instability would then 
fall in the RII region. 

Additional information could come from directly detecting the 
accretion disc itself in quiescence. A trapped disc will have both a 
higher surface density and temperature than a standard accreting 
disc with the same rj n . For EX Lupi in quiescence (using the pa- 
rameters given above), our model predicts a surface temperature 
from internal viscous dissipation of 530-650K (depending on m), 



compared with 290-570K for a Shakura-Sunyaev disc. This mod- 
est temperature increase is likely insignificant compared to heating 
from the central star, which will thus stil l determine the gas be - 
haviour at these low accretion rates (e.g. [Muzerolle et alj|2004r) . 
The surface density in the disc will also increase modestly, from a 
maximum of between 3-40 gem" 2 (p ~ 4 x 10 -11 - 6 x 10~ 10 
gem" 3 to 40-80 gem" 2 (p ~ 4 x 10~ 10 - 10~ 9 gem -3 ). 
At these low densities the disc is probably optically thin (e.g. 
iDullemond & Monnien l20l(il) . so that it is likely not possible to di- 
rectly distinguish between a trapped disc and normal accreting one 
except through detailed modelling of spectral features in the disc. 

Unlike the accretion instability (which only occurs under cer- 
tain conditions in the disc-field interaction), the build-up of mass 
just outside the co-rotation radius should be generic in stars with 
strong magnetic fields and relatively low accretion rates. The best 
opportunity to detect a trapped disc directly would come from a star 
with a strong dipolar magnetic field, high spin rate and low accre- 
tion rate, so that the disc would be optically thick rather than thin 
(as would be expected for the given accretion rate), and could be 
directly detected through modelling of its spectral features. 

5.2 NS transients w/ weak recurrent outbursts 

In DS11 we suggested that the large ratio between the spin-down 
and viscous timescales in accreting millisecond pulsars makes it 
unlikely that the disc will become trapped, since as rh drops r- m 
moves outwards much faster than r c . This agrees with the typical 
outbursts seen in transient X-ray pulsars, which show outbursts of 
variable duration followed by long periods of quiescence and are 
believed to b e triggered by ionization instabilities in the disc (e.g. 
lLasotall200lh . 

However at least tw o accr eting X-ray pulsars, NGC 
6440 X-2 dHeinke et alj l2010h and IGR J00291+5934 
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dHartman. Galloway & Chakrabartvl 1201(f) have been observed 
to undergo weak, recurrent outbursts on much shorter timescales 
(around 30 days) than would b e predicted from an ionization- 
instability model dLasotd I200TI) . The weakness of the outburst 
has been interpreted to mean that the disc is never completely 
ionize d during the outburst lHartman. Galloway & Chakrabartvl 
(2010), which suggests the ratio of Tsd/TUsc becomes small 
enough to trap the disc as when the accretion rate drops. This large 
reservoir of matter would then be present when the accretion rate 
again increases, and so could allow more frequent outbursts. These 
sources thus offer an excellent opportunity for comparison with 
our model of a trapped disc. 

IGR J00291+5934 is a 599Hz X-ray pulsar in a 147 minute or- 
bit, which historically has shown standard single outbursts (in 1998, 
2001, and 2004) followed by long periods of quiescence. How- 
ever, in 2008 it underwent a substantially different outburst, with 
two weaker outbursts (respecti vely lasting nine and 13 days) sepa- 
rated by 27 days of quiescence ( lHartman. Galloway & Chakrabartvl 
l201fj|) . The authors estimated that the total mass flux of weaker 
outbursts together was approximately equal to the single 2004 out- 
burst. Similar outbursts have been seen in NGC 6440 X-2, a re- 
cently discov ered X-ray binary with a frequency of 442Hz in a 57 
minute orbit dAltamirano et alj|2010h . Since discovery on July 28, 
2009, five distinct outbursts have been observed, each typically last- 
ing around four days at a luminosity of around Lx ~ 10 35 erg s _1 , 
with a minimum recurrence time of around 31 days. 

The short recurrence timescales suggest that a considerable 
amount of mass is stored in the disc when the accretion rate drops, 
which can then fuel another outburst when rh increases again. In 
both sources the quiescent X-ray luminosity is at least 4 orders of 
magnitude lower than the outburst luminosity, indicating a large 
change in accretion rate (although the X-ray luminosity may not 
be a good indicator of bolometric luminosity at such low accretion 
rates, so there is some uncertainty in the change in rh). 

During two outbursts, the power spectrum of NGC 6440 

X-2 was dominated by a strong low-frequency QPO ( ~ 1 Hz) 

dPatruno et al]|2010l: lHartman. Galloway & Chakrabartvll20l61) . A 
similar QPO was previsouly d etected in the tail an outburst in X- 
ray pulsar SAX J1808.4-3658 dPatruno et ai]|2009l) which was at- 
tribute d to the disc instability presented in DS 10 and lSpruit & Taarr] 
( 1993). A much weaker QPO at 0.5 Hz was observed in the out- 
bursts IGR J00291+5934 and might also be attribut ed to the disc 
instability (Hartman, Galloway & Chakrabarty 2010). 

The fact that all of these QPOs were observed at a similar 
frequency presents a challenge to our interpretation. The outbursts 
seen in SAX J1808. 8-3658 follow the more conventional pattern 
of a large outburst followed by a long quiescence, which would 
suggest that a dead disc does not form. However, the similarity of 
the QPO frequency (which is roughly correlated with the viscous 
timescale near n n in our model) suggests that the viscous timescale 
is similar in all these sources. This contrasts with our suggestion 
that the viscosity in these weak outbursts is substantially smaller 
than in larger outbursts. The number of uncertain parameters in 
our model (in particular this assumes that Ar2/Vi n is similar in 
all sources) means that the picture we have presented could still 
be valid, but then the similarity in QPO frequency in both (or all 
three) of these sources is probably coincidental. On the other hand, 
a strong low-frequenc y QPO has currently be en observed in only 
two sources, while as lAltamiran o et all d20ld) have noted, weakly 
outbursting sources like NGC 6440 X-2 will generally be missed in 
current surveys. 

In our present model, the disc-field interaction is not 



responsible for directly regulating the accretion rate onto 
the star to control the outbursts, as was suggested by 
lHartman. Galloway & Chakrabartvl d20ld) . The field can temporar- 
ily halt accretion onto the star, but only on the much shorter 
timescale on which the 1Hz QPO is seen. Thus the mechanism that 
directly triggers an accretion outburst remains undetermined. 

5.3 Persistent X-ray pulsars 

A small subset of accreting pulsars with low-mass companions 
have sufficiently persistent X-ray emission to measure the spin 
change directly and study the connection between spin change 
and luminosity. Three well-studied persistent sources (HerX-1, 4U 
1626-67, and GX 1+4) all show episodes of spin-up and spin-down, 
althou gh the pattern of sp in-up and spin-down is different in all 
three dBildsten eta"i1ll9970 . A fourth, (GRO J1744-28) shows clear 
spin-up when in outburst, allowing the relationship between v and 
rh to be constrained. 

The binary properties of all these sources varies consider- 
ably, both in their companions and orbital periods. However, it is 
still possible to make some general observations that apply to all 
sources. First, there is a clear correlation between luminosity and 
spin change, with higher luminosity corresponding to spin up, and 
vice versa. Second, when the star is spinning up, the spin-up rate 
is substantially lower than would be predicted from the X-ray lu- 
minosity. Third, where the relationship between v and rh can be 
measured, v oc m 7 , where 7 > 6/7, i.e. the scaling expected 
for n n oc ra~ 2//7 . Finally, the spin-up torque is often compara- 
ble (within a factor 2) to the spin-down torque, with only a change 
in sign. 

The first observation is universally expected in magneto- 
spheric accretion, since in order to spin up the star the matter must 
be able to overcome the centrifugal barrier imposed at r c . The sec- 
ond point is also readily understandable in our model, which as- 
sumes that there is a transition region for n n ~ r c in which there is 
both accretion and spin-down torque. This behavi our is also seen in 
MHD simulations of the disc-star interaction (e.g. lRomanova et"al] 
120041) . For high accretion rates, our model predicts v oc rh 1 , 
which is close to the dependence measured in some X-ray pul- 
sars in outburst. For example, the 1995/1996 outburst of GRO 
J 1744-28 (a disc-fed binary), was measured to follow i> oc rh ' 96 
(Bild sten et "all ll997l) . This is closer to observation than the stan- 
dard model, which predicts an index of 6/7. 

The final point is most puzzling in a model in which ri n oc 
rh~ 2 ' 7 , because it requires that the change in accretion rate is finely 
tuned to preserve this symmetry. It is more plausible for our model, 
however, which shows a steep decrease in i> around rh/rh c = 1, 
and then a flattening around the minimum in v. Thus, a small de- 
crease in rh will cause an initial steep transition to spin-down be- 
fore settling around the minimum in v. There is no comparable 
turnover at high rh, but in order to switch between spin-up and 
spin-down, the star must have an average accretion rate such that 
»"in ~ r c , while in order to be a persistent source, the accretion rate 
at large distances cannot be very variable. If we assume that the av- 
erage accretion rate has been steady enough on long timescales to 
affect the star's rotation rate, then there will be a natural spin-down 
torque, set by the minimum in [2] and we would expect that on av- 
erage the accretion rate would fluctuate around the level that keeps 
ri n close to r c . This still does not account for the observations (for 
example of 4U 1626-27) in which the magnitude of the torque is 
very nearly symmetric. 

Our present model is too simple to make quantitative com- 
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parisons to data. For one thing, the relation between v and m 
depends on the least certain parts of our model: both the tran- 
sition widths Ar/nn and Ar2/n n and the connecting functions 
y m and y^. There is also much more variation in the luminos- 
ity at which the system moves from spin-up to spin-down than 
is generically predicted in any simple magnetospheric accretion 
model. This variability could indicate fluctuations in (rf), and 
Ar/ri n, which would both affect th e magnitude of spin-down 
torque. Agapitou & Papaloizou (2000) demonstrated that the value 
of Btf, changes if field lines are able to diffuse radially through the 
disc. The d iffusion of field lines is expected to happen on viscous 
timescales dFromang~& Stone 2009), so that the details of the disc- 
field interaction could in fact rely on the complicated interplay be- 
tween the accretion rate and the field strength and geometry. Such 
questions require detailed MHD modelling which are beyond the 
scope of the present work. 



6 DISCUSSION 

At low accretion rates, the strong magnetic field of a star can 
strongly alter the structure of an accretion disc, preventing it from 
either accreting or expelling material outward, and preventing the 
inner edge of the disc from moving considerably away from r c . As 
we have discussed in this paper, this disc state can lead to substan- 
tially different behaviour from what is conventionally assumed, and 
could explain a variety of observed phenomena in magnetically ac- 
creting stars. However, in order to construct our model we used a 
parameterized description for the interaction between the magnetic 
field, introducing the free parameters Ar/n n , and Ar2/n n . Here 
we discuss what sets the value of these parameters, and whether 
these parameters could vary in time, which could account for some 
of the additional behaviour not described by our current model. 

The first parameter, Ar/r- m is the width of the direct interac- 
tion between the field and the disc in the disc's innermost region. 
We have assumed in all cases that Ar/nn < 1, so that the inter- 
action region (where the field lines remain closed long enough to 
exert a significant stress on the disc) remains fixed. MHD simu- 
lations suggest that this region will fluctuate in time as field lines 
open and reconnect, but we have assumed that this will occur on 
timescales of ~ P», which are much shorter than the viscous 
timescale of the instability (T v i sc > 10 3 P„). Numerical MHD sim- 
ulations do sometimes show variability on longer timescales (e.g. 
Rom ano va et al. I J2009k which could lead to changes in Ar/nn 
between individual outbursts. Given that even small changes in 
Ar/nn can significantly change the appearance and amplitude of 
the outburst (seen in figure [6), variability in the disc-field interac- 
tion could produce variability in the appearance of the instability, 
especially if the instability occurs in RII. 

A second possible source of variability could come from the 
radial evolution of the magne tic field itself. This was studied by 
( Agap itou & Papa loizou 2000), who found that radial field-line dif- 
fusion through the disc can reduce the effective torque at n n , and 
increase the extent of the closed-field line region. The timescale 
for magnetic diffusion is typically assumed to be about the same 
as for viscous diffusion (whic h is supported by recent MRI calcua- 
tions; Fromang & Stone 2009), so that the radial diffusion of closed 
field-lines outward could compete with accretion of matter (and 
field lines) inward, creating an additional, longer timescale source 
of variability in Ar/n n . This could allow the instabil ity to appear 
and disappear at the same m (as has been seen in IPatruno et alj 



2009, where the 1-Hz QPO was seen in portions of the outburst but 
not others at the same flux). 

The physics that sets Ar2 /n n (the transition length between 
the accreting and non-accreting solutions) is considerably less 
clear, although presumably the same variations in the disc-field in- 
teraction that change Ar/n n could also change Ar2 /n n . In DS 1 1 
we discussed the model proposed bv lPerna. Bozzo & Stella|j2006l) , 
which considers a strongly inclined dipole, so that the disc trunca- 
tion radius lies both inside and outside r c , allowing for simultane- 
ous accretion and confinement of material. If this is the case then 
Ar2 / n n should remain constant for a given system, and the vari- 
ability must come from fluctuations in Ar/n n . 

Another open question remains what regulates the mean ac- 
cretion rate onto the star. In all our simulations we have treated the 
accretion rate, m as a free parameter, and shown how the instability 
and spin-evolution of the star change as a result. However, obser- 
vations of neutron star binaries suggest that the accretion rate from 
the donor star remains roughly constant Bi ldsten et alj d 19971) . and 
variatio ns in m are o ft en assumed to arise from ionization insta- 
bilities jLasota|[200lh . I Hartman. Galloway & Chakrabartvl d201Ch 
suggested that the disc-field interaction could regulate the accre- 
tion rate onto the star by storing it up, in a similar way to what 
we have discussed here. However, our model does not agree with 
this conclusion: the disc-field interaction can only halt accretion 
long enough to build up mass to drive the disc instability and only 
substantially alters the inner regions of the disc. Our trapped disc 
solutions evolve as a response to the decreased accretion rate in the 
disc, and only become relevant once the average accretion rate has 
decreased below a threshold, m CI i t /rh c ~ 4Ar/n n (sec. 13. It . 

On the other hand, if radiation feedback from accretion onto 
the star plays an important role in regulating the viscosity (and so 
accretion rate) in the disc, then temporarily suppressing accretion 
onto the star might cause the level of ionization to drop and drive 
the source back into quiescence. 

A final question remains: why do some sources show the ac- 
cretion instability and trapped disc behaviour and not others? The 
vast majority of outbursting neutron stars show long periods of qui- 
escence, consistent with being driven by the ionization instability 
with ri n far from r c , rather than remaining trapped near r c . Of 
the two sources we have considered, IGR J00291+5934 has pre- 
viously shown a normal outbur st pattern (a much stronger outburst 
follow ed by a long quiescence: lHartman. Galloway & Chakrabartvl 
l20ld 

In their discovery paper for NGC 6440 X-2, Alta mirano et alj 

d2010h remark that the short recurrence time and weak out- 
burst likely cause similar NS sources to be missed by cur- 
rent surveys, so it is possible that such sources are common 
but remain undiscovered. As well, in DS11 we showed that 
the very small ratio T v i ac /TsD ~ 10 -17 for rapidly spinning 
pulsars made it more likely that the disc would become un- 
trapped (with n n moving rapidly away from r c ) as rh — > 0. 
lHartman. Galloway & Chakrabartvl d2010h suggested that the disc 
might remain weakly ionized during the entire outburst, which 
would then decrease the viscosity and increase T v i sc , perhaps 
enough to prevent the disc from becoming untrapped. If the RI in- 
stability is responsible for the 1Hz QPO seen in NGC 6440 X-2 
and SAX J1808. 4-3658, then the fact that only sometimes mani- 
fests is consistent with the idea that Ar/n n is variable (since the 
instability only occurs over a small range in Ar/n n ). 

For T Tauri stars, the main constraint on the instability is the 
presence of a strong (~ 1 kG) dipolar magnetic field. Although 
most T Tauri stars show a strong surface field component, the dipo- 
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lar component can vary considerably, to the point that it is not clear 
whether the disc-field interaction considered here is primarily re- 
sponsible for the low spin-rate of the star. The inclination between 
the rotation axis of the star and the magnetic field could also af- 
fect the instability by changing A^/V™ (as above). This would 
suggest that, as a class, EXors are T Tauri stars with moderate ac- 
cretion rates and strong, aligned dipolar magnetic fields compared 
with normal T Tauri stars. 



7 CONCLUSIONS 

As found in DS10, accretion from a disc onto a magnetosphere 
tends to take place in what we have called 'trapped' states. The 
inner edge of the disc hovers around the corotation radius in this 
state, even at very low accretion rates. The accretion in this state 
can be steady or cyclic. We have investigated here the conditions 
under which the cyclic form takes place, the kinds of cycles, their 
amplitudes, cycle periods and their effect on the average torque on 
the accreting star. 

Two forms of cycle are found; the first form (labeled 
RI here) is the accretion/ dead disc cycle already suggested by 
Sunyaev & Shakun] d 1977h . A pile-up of mass at the 'centrifugal 
barrier' outside corotation is followed by an episode of accretion 
emptying the pile and the beginning of a new cycle. Mass loss or 
'propellering', whether it takes place or not, is a separate piece of 
physics not needed for understanding the effect of a centrifugal bar- 
rier on a viscous disc accreting on a magnetosphere. 

Whether accretion is cyclic or continuous is found to depend 
on the width of the transition from an accreting state (inner edge in- 
side corotation) to a suspended accretion state (edge outside coro- 
tation). A narrow width leads to RI cycles, a broader transition 
to continuous accretion. The effect of these cycles is to increase 
the average spindown torque on the star. A second form of cycle 
takes places when the average accretion rate is close to the stan- 
dard accreting case (edge inside corotation). It decreases the aver- 
age torque on the star. 

A review of the available literature shows forms of variability 
that may be related the the cyclic forms of accretion found here, 
the most promising ones may be the 'EX Lup' outbursts, and a new 
form of cyclic accretion found in X-ray binaries NGC 6440 X-2 
and SAX J1808.4-3658. 
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